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ABSTRACT 


Optimum resolution was determined for a series of 
straight chain hydrocarbons under time normalization. 
The variable parameters were programmed temperature rate, 
isothermal temperature, normalized net retention time, 
hes Diow rate, and: column length, The percent liquid 
phase in each column was such that the total amount of 
liquid phase remained constant for each column in a series 
or varying column lengths. The results showed that if 
one worked with an amount of liquid phase normally packed 
into a one meter column using 25 percent liquid phase 
packing, the resolution varied by ieéss than one unit 
when the column length varied from approximately 1.5 to 
5.O meters for the separation of a given pair of straight 
Chain hydrocarbponss ‘lor a giver helium Tlow rate and for 
a given programmed temperature rate. For column lengths 
lees than approximately 1.5 meters the resolution decreased 
Sharply wich decreasing column Jeneth.. The optimum pro- 
grammed temperature rate was dependent on the compounds 
being examined e.g. as the straight chain hydrocarbon 
number increased the optimum programmed temperature rate 
decreased, For all the normalized net retention times, 
column lengths and hydrocarbon numbers, the resolution 


increased with helium flow rate. 
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PART Ht 


INTRODUCTION 


In time normalization the parameters of the gas 
chromatographic system are such that the last component 
of the separated sample comes” OUT’ Or “the column. With” 4 
fixed net retention time as one or more parameters are 
Changed. This gas chromatographic technique has received 
relatively little attention in the field of gas-liquid 
chromatography. The first experimental work using time 
normalization was done by Karger and Cooke (1) who showed 
that for maximum resolution under normalized time conditions 
an optimum column length exists for both packed and 
capillary columns. The column used was a PAO ie Ord wy 
copper tubing packed with 30/60 mesh Chromosorb PH coated 
With S percent. Dow’ 710 silicone oil.” The sample Injected 
was oe of 1 percent (1 to 1) dodecane and tetra- 
decane in n-hexane. 

The optimum length for maximum resolution for 
normalized retention times of 1200 * 60 seconds and 
2860 + 90 seconds for tetradecane was 30 feet and 55 
feet. The optimum length for maximum resolution under 
isothermal conditions of 130°C was 80 RSE 

In the same experiment a capillary column was pre- 


pared by evaporation of a 10 percent solution of silicone 
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Oil. When both capillary and packed columns were norma- 
lized to a retention time of 1200 seconds the optimum 
length for maximum resolution were 150 feet and 30 feet. 
In another publication (2) by the same authors the 
particle size and average carrier gas velocity were also 
optimized under normalized time conditions. The effect 
of mesh size on resolution was examined under constant 
inlet pressure and constant average carrier gas velocity 
under time normalization. Five columns of Chromosorb P 
with mesh sizes 45/60, 60/80, 80/100, 100/120, and 120/140 
each coated with 4 percent Dow 710 silicone oil were 
prepared. When a solute such as 2-methylnaphthalene 
with a low relative retention dependence on temperature 
or a solute such as n-octane with a high relative retention 
dependence on temperature was injected it was found that 
large particles (45/60 mesh) were best for low inlet 
pressures (low oon rates) and small particles (120/140 
mesh) were best for high inlet pressures (high flow rates). 
In a study where average velocity was kept constant 
they used a column of constant length and 4 percent Dow 
{10 silicone oil. When solutes such as the methyl- 
naphthalenes were injected under isothermal and nor- 
malized time conditions, resolution increased with de- 
creasing particle size until 80/100 mesh. With further 


decrease in particle size the resolution remained constant. 
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When solutes such as the iso-octanes with a high relative 
retention dependence on temperature were injected under 
constant average velocity and normalized time conditions 
resolution was about the same for large (45/60 mesh) and 
for small (120/140 mesh) particles. 

For normalized times of 134 and 186 seconds, the 
resolution for mesh sizes 45/60 and 120/140 were almost 
equivalent, while for the normalized time of 94 seconds 
the resolution for mesh size 45/60 was 9 percent higher 
than for mesh size 120/140. The authors concluded that 
under constant average velocity and normalized time con- 
ditions "a better resolution is possible without loss in 
retention time by operation above the velocity for Hae 
with a simultaneous reduction in temperature" (2). 

Grushka, Yepes-Baraya and Cooke (3) showed the effect 
of temperature, carrier gas velocity and amount of sta- 
tionary, phase on resolution under normalized .time con- 
eaciee Packed columns used were 2m x 1/8 in., o.d., 
copper tubing, packed with 3 to 30 percent w/w Apiezon L 
on 80/100 mesh AWDMCS Chromosorb W. Capillary columns 
were made from a 316 stainless steel 0.01 in., i.d., 150 
feet in length and coated with 17 percent w/w Apiezon L 
in toluene. When experimental runs with this length of 
column were completed, the column was cut in such a way 
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sample mixtures used for the packed columns were: 
a) 1 ce of 3-methylnonane and 1 ce of 5-methylnonane 
diluted to 25 ce with acetone; b) 1 ce of sec-butylbenzene 
and 1 ce of tert-butylbenzene diluted to 25 ce with benzene. 
Mixtures used for the capillary columns were 3-methyl- 
nonane and decane and a mixture of l-methylnaphthalene 
and 2-methylnaphthalene. iis amounts of these solutions 
were injected. 

micreasinge “percent Tidquid pnase at’ constant temperature 
and velocity for the 3-methylnonane/decane mixture on 50 ft 
Capillary “columns showed that an optimum range of liquid 
Prin “circ kness dp exis vo ln "wile the resolevlon rs maxki= 
mized (d, supa Pncreasing percent’ Frqwtd'phaseon 
packed columns at constant velocity and temperature showed 
that the increase in resolution with percent liquid phase 
Wae novU-~as "great -as or *thecapi liary “columns >*-The reso- 
lution of the methylnonanes was approximately the same 
on all their columns while the resolution of the butyl- 
benzenes showed a small maximum. 

Under time normalization for the decane/3-methyl- 
nonane mixture on capillary columns at constant velocity 
and. retention time of 211.3 seconds the results showed 
that the resolution decreased with increasing percent 
liquid phase and temperature in spite of an initial de- 


crease in H (minimum at d, = Tae for decane. For packed 
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columns under time normalization for the 3- and 5-methyl- 
nonanes the results showed that the resolution also de- 
creaséduwithsineneasingspercenthliquidecphases «Theyncon- 
cluded that resolution is best for the lower loaded column 
operated at a lower temperature. 

Under time normalization where velocity instead of 
temperature is changed for both capillary and packed 
columns the results showed that an optimum range exists 
in the amount of stationary phase and in the velocities 
of the gas phase. For capillary columns the optimum 
film thickness was approximately ea with the carrier 
velocity about 47 cm/sec. For packed columns maximum 
resolution occurred at about 5 percent loading at about 
7.9 cm/sec. Their theoretical equation indicated a 
maximum resolution at a liquid film thickness of 0.3 
and a carrier velocity of about 46.5 cm/sec for the 
capillary column and a maximum resolution a 3.5 percent 
liquid loading and carrier velocity at about 6 cm/sec 
for the packed column. In all the cases above the optimum 
liquid phase loadings and velocities were simultaneously 
at the optimum. 

Grushka, Yepes-Baraya and Cooke (3) derived an 
equation from R = 1/4(L/H,) 7? (a-1/a) (ko/Ko+1) in which 
they substituted a simplified van Deemter equation for 


2 a ‘ 
Hii 2yD,/u os qdek,u/D, (k,+1) where L = column length, 
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frgpet impos = 2 seadt i i ielle 


Ome relative metymetentilon wime, Ky ey capaci by) ratdo-vior 
umeneeconu peak, y= a. correction LTactor to allow for 

the sinusoidal paths traversed by the carrier gas in the 
column (also known as the tortuosity factor or obstruction 
factor and thaving a’ value of 0.5 - 0.7), D, = the dif- 
fusion icoetfi cient tof tthe: isolipe “in the gas phase, u = 
the carrier gas average velocity, q =a geometric factor 
known as a configuration constant the value of which 
depends upon the shape of the liquid phase pool (assuming 
a uniform film thickness, q has a value of 2/3), d, = the 
film thickness, and Dd, =the -dirrusion coefficient of the 
solurce in the liquid phase. 

The assumptions made were: a) the eddy diffusion 
term and the resistance to mass transfer i.e. the radial 
(gas resistance) diffusion term are vanishingly small; 

b) the capacity ratio is proportional to the percent 

liquid phase, i.e. k, = a(%); and c) the film thickness 

d. is given by d, = 2/S 2, where % = the percent of liquid 

phase, on the packing, S = the surface area per gram of 

the support, and Jae the density of the stationary phase. 
substituting the simplified van Deemter equation 

for H, into the equation for R and then substituting the 
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equation for u = ju, = L(1+k,)/t (where tao = the nor- 


Re 
malized retention time for the second peak, u = the average 


Fas VeELloci cy, a the gas velocity at the column outlet 
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and j = the compressibility correction factor) one then 


obtains 


WAIL (k5/k5+1) (a-1/a) 


es 
2yD, tro i 2k,L(%) 


ot) SD eer Ss A) 


which reduces to 


1/4 {L k,(a-1/a) 


2yDatpo(Kot1) ' 2k, L(K +1) (%) © 


2 
A 3D) to (SK) 


Using k, = a(%), R becomes 


1/4 [L (a%) (a-1/a) 


2D tao (age) . oLa(%)>(a%+1) 
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E 3PytRo(S4) 


which is the equation derived by Grushka, Yepes-Baraya 
and Cooke (3) and was used to describe experimental 
results for packed columns under normalized time conditions. 

The columns used in previous experimental work were 
either of constant length (2,3) or where various column 
lengths were used (1) the amount of liquid phase inside 
the column increased with column length. In my work 


attention was focused on a series of columns of varying 
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lengths, each column containing the same amount of liquid 
phase. Time~normalization was applied to the determination 
of maximum resolution when column length, helium flow 

rate, isothermal temperature, and programmed temperature 
rates were varied for’a sample containing a series of 


straight chain hydrocarbons. 
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PART 1 
EXPERIMENTAL AND EQUIPMENT USED 
1 EXPERIMENTAL 

The liquid phase selected was UCW98. The solid 
support was unsilanized Chromosorb WAW 60/80 mesh. The 
column was 1/8 in., 0o.d., Anaconda soft copper tubing. 

A 5%, 8.3%, 25%, and 50% liquid loading on solid 
support was prepared by dissolving 0.608, 1.00g, 2.96g, 
and 6.02g UCW98 in diethyl ether and added to 12.07g (50.29 
m1), 12.03g (50.14 m1), 11.85¢ (49.37 m1), and 12.05¢ (50.21 
m1) of Chromosorb WAW 60/80 mesh. Four columns of lengths 
5.0m, 3.0m, 1.0m, and 0.5m were cleaned with acetone, al- 
lowed to dry and were filled with the 5.0%, 8.3%, 25.0%, 
and 50.0% packing using the gravity method. In this 
method the dry coated support is slowly poured through 
a funnel into the vertical tubing while gently tapping 
or vibrating the column (4). The columns were then bent 
into 8 inch diameter coils. The 5.0m column was made 
from two 2.5m lengths and were connected by 1/8 inch 
Swagelok unions. Quartz wool was used to plug the ends 
of the copper tubing. Each column then contained 0.144¢ 
of liquid phase, calculated as (vol. of column / vol. of 
packing prepared) x wt of liquid phase used to prepare 
the packing. An equivalent set of four columns was pre- 


pared and used as reference columns. Before use each 
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column was conditioned at 200°C for 10 hours during 
which a helium flow was maintained. After conditioning, 
approximately 40 mw of a GLC Column Conditioner (Silyl1°8, 
Pierce Chemical Co... P.O, Box 117, Rockford, Illinois.) 
was injected into each column at 200°C. The treatment 
was repeated once after a 5 minute interval. 

A ae Hamilton syringe was used to inject a sample 
size of ae ee n-nonane, n-decane, n-hendecane, and n- 
dodecane with relative peak areas 1.0, 1.57, 2.56, and 
4.59. The same sample size with the four components 
mentioned above was used throughout the experiment. The 
gas chromatographic peaks were digitized using a digital 
integrator with printer. The digitized peaks were then 
analyzed by an IBM 360/67 computer which then calculated 
HETP and resolution using a Fortran program which appears 
at the end of this report. 

Helium flow rates were measured with a soap bubble 
flow meter at 25°C and were corrected for water vapor 
pressure in the flow meter by the formula 
V. = V_(P_-P Wace where V, = corrected flow rate at 
outlet column pressure, Nhe = flow rate measured with 
soap bubble flow meter at outlet column pressure, E. a 
atmospheric pressure, and Pa = the vapor pressure of 
water at 25°C. All retention times were measured with 
outlet column at atmospheric pressure and were not 


corrected with the compressibility factor Jj. 
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Time normalization, coupled with a wide variation 
of programmed temperature rates and a wide variation of 
net retention times necessitates a wide variation of 
initial temperatures for the column. However with the 
temperature programmer used in this study, initial temp- 
eratures were restricted to those above room temperature. 
The large size of the oven used restricted the heating 
rate to less than 14 deg/min. These restrictive conditions 


resulted in some graphs which appear to be incomplete, 


2 EQUIPMENT USED 
The gas chromatographic system was constructed from 

separate units-ot equipment... A schematic diagram of the 

system isigivem_in Figure 1. 

AL ebrooks i Llowicontrollers:. 

B. Injection block - one in which the helium passed 
close to the septum. 

C.. He Kul». pWOmMOde—ControlLlem el AePoi-rinevrunenge: Wen 
It controlled the temperature of the injection block 
au 2007 C. 

D. Programmer - Varian Aerograph Linear Temperature 
Programmer Model 330. 

Heo Commmn or 1/6 Ing, (0.0, ,,. copper. tubing or varying 
lengths packed with UCW98 liquid phase (Hewlett- 
Packard) on Chromosorb WAW 60/80 mesh (Hewlett- 


Packard, Avondale Division). 
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Fig. 1 The gas chromatographic system 
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Gow-Mac power supply, Model 40-05-C, and Gow-Mac 
thermal conductivity cell, Model 9285-D, Gow-Mac 
Liscrumeneun CoLGwMadison, New Jersey. The detector 
current was 180 ma. 

Aerograph digital integrator - Model 471-42 - manu- 
facwuredviby-tintotvronics? Comprs§ Houstvorg? Texasy 
Victor digit-matic: printer - Victor Comptometer 
Corp: ,UChicaso Cea nois . 

Filter - Trap number 236 - Guild Corporation, Bethel 
Par ig,ueraner Ua, A 

Temperature Controller - RFL Model 70-115/208/230V - 


RUuPsiAY Industries ,q@.lnc .,aBoontonienNew Jersey mn Uses A. 
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PART 21 
RESULTS AND DISCUSSION 


i” HETP vs Helium Flow Rate 

Fig. 2 shows the effect of isothermal temperature, 
column length, straight chain hydrocarbon number, and 
helium flow rate on HETP, The three isothermal tempera- 
tures selected were 05°C, 110°C andr 1265°C: 7 in alithe 
cases shown the HETP decreased with isothermal temperature 
from 125°C to 95°C. For a given isothermal temperature 
increase, the percent increase in HETP increased as the 
column length increased and as the Straight chain hydro- 
carbon number decreased e.g. for a change of 15°C in the 
isothermal temperature, the change in HETP for n-nonane 
was dramatic while for n-dodecane the change in HETP was 
insignificant. This can be explained as follows. HETP, 
column length, resolution, the relative volatility a, and 
the capacity factor k are related by the equation (1) 
HETP = (L/16R*) (a-1/a)°(k/«+1)? ; “he-ratios:d=-7/a ard 
k/k+1 change most when a is close to 1 and when k is 
small i.e. when the column temperature is near the boiling 
point of a sample component. [In this experiment the 
column temperature was nearest the boiling point Hf n- 
nonane and farthest from the boiling point of n-dodecane, 
For a given column length e.g. 5.0m and e.g. isothermal 


temperature at 125°C and helium flow rate 21 ml/min, 
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Pig. 2 HETP vs Helium flow rate. Columns from left 
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Vora ale On. COLUM Tene wie OL UO, oO, ule. ves eile 
5.0 meters. Rows from top to bottom are for n-nonane, 
n-decane, n-hendecane, and n-dodecane. A, B and C are 
PoOuNelMal vemperalvures Git leo. Cw iO. UC andwO>e Ce 

the HETP decreased as the hydrocarbon number increased. 
nee resoOluL tone 8b 125 °C tor 


C C and s0n 2 —C 


Os aCe ae brell 1 ane 
Werens 106.14). 06-and 502. othe capacity ratios at 125°C 
for Cg» Cro» Cyz, and Cj, were 0.7, DP Ouete ¢ oom: And wD 
The ratios k/k+1 at the same temperature for the same 
series were 0.411, 0.561, 0.698, and 0.808. The changes 
Ak Ake 1) for Cg-Cr9; Cyo7Cq, and C,,-C,5 were 0.150, 0.137 


and 0.110... The change in resolution AR between C.-C 


OF a0 
end (Ca n=O was 1.60 and between Cc CG BC Ct. ae 


LOisda: Oty 1h ial t 1e 
was 0.16. The value of a for all the hydrocarbons was 
approximately 1.8 Since temperature was constant and the 
change Aa from one hydrocarbon to the next will be con- 
Sidered negligible. Therefore the change in HETP will 
increase as. the hydrocarbon number. decreases, i.e. the 
HETP changes most for n-nonane and least for n-dodecane. 
This was experimentally observed in Figure 2, 

With a few exceptions the absolute value of HETP 
increased as the column length increased from 0.5 meters 
to 5.0 meters. For the exceptions, which were 95°C for 
n-nonane, 95°C and 110°C for n-decane and 95°C, 110°C 
and 125°C for n-hendecane, the HETP decreased slightly 


(less than 0.3mm) when the column length increased from 


O.5 meters to 1.0 meter and then the HETP increased signi- 
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ficantly when the column length increased further from 
1.0 meter to 5.0 meters, 

The optimum helium flow rate increased with column 
length and was independent of the straight chain hydro- 
carbon number. For example the optimum helium flow rate 
for the’ O.5 meter column for. m-nonane; n-decane, n=-hendecane, 
and n-dodecane was approximately 13 ml/min and for the 
5.0 meter column for the same straight chain hydrocarbons 
the optimum helium flow rate was approximately 30 ml/min. 
As the straight chain hydrocarbon number increased the 
optimum flow rate became less definable. The HETP plot 
for n-dodecane for example was almost independent of 
helium.flow rate. 

From the van Deemter plots in Fig. 2 one can deduce 
that the optimum parameters are a column length of approx- 
imately 1 meter, an isothermal temperature of 95°C and 


a helium flow rate of approximately 20 ml/min. 


2 Resolution vs Helium Flow Rate 

Fig. 3. shows ‘the effeet® of? helium Flow rate, column 
length and isothermal temperature on the separation of 
n-nonane, n-decane ,, n-hendecane... and, n=dodecane.). In tall 
the cases shown the resolution increased as the tempera- 
ture decreased from 125°C to 95°C at 15°C intervals. 
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Pig. wes REeSOluLion.vs helium flow rate... Columns from 
Left to rvent are Lor column -fereurs oreo ys Piao ees 
and 5 .Osemeters... Rows .trom,.top.to bottom are for 
n-nonane-decane, n-decane-hendecane and n-hendecane- 
Gogecane.. Dow andy are Lor lecothermal téemperavures 
Or 95 €y 110°C and. Ie5 
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percent increase in the resolution increased as the 
columm length increased andsas the straight-chain hydro- 
carbon number decreased e.g. fora «change of)15°C in the 
isothermal temperature, the change in resolution for 
n-nonane-decane for a given column length was greater 
than the change in resolution for n-hendecane-dodecane 

at the same column length. For all three isothermal 
temperatures and all the straight chain hydrocarbons the 
absolute value of the resolution increased with column 
length from 0.5 meters to 5.0 meters although there was 
no significant increase when column length increased from 
3.0 meters to 5.0 meters where it appears to have reached 
a limiting or optimum value. 

The optimum helium flow rate increased with column 
length and was independent of the straight chain hydro- 
carbon number. For example the optimum helium flow rate 
for the 0.5 meter column for the separation of n-nonane- 
decane, n-decane-hendecane and n-hendecane-dodecane was 
approximately 13 ml/min and for the 5.0 meter column for 
the same straight chain hydrocarbons the optimum helium 
flow rate was approximately 30 ml/min. 

For the separation of a given straight chain hydro- 
carbon, i for a given.column length and: for a given. iso- 
thermal temperature the resolution changed by less than 
one unit when the helium flow rate varied from approxi- 


mately 10 ml/min to 36 ml/min. 
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From the plots in Fig. 3 the optimum parameters are 
a column length of about 4 meters, an isothermal tempera- 
ture of 95°C and a helium flow rate of approximately 
20 ml/min. The optimum temperature and optimum helium 
flow rate are in agreement with the van Deemter plots in 
Fig. 2... However the resolution :vs helium flow rate.plots 
in Fig. 3 indicate an optimum column length of 4 meters 
while the van Deemter plots in Fig. 2 indicate an optimum 
column length of 1 meter. Here we have a case where 


HETP and resolution both increase with column length. 
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3 Time Normalization With Net Retention Time lor 
N-Dodecane Fixed at 300 t 5 Seconds 

Figs. 4, 5 and 6 show the effect of helium flow rate, 
column length, isothermal temperature, and programmed 
temperature rate on the resolution of n-nonane, n-decane, 
n-hendecane, and n-dodecane when the net retention time 
for n-dodecane was eee ae 300 5s 5 seconds. 

The three helium flow rates used were 6.8, 13.5 
and 20.3 ml/min. In all the cases shown the resolution 
increased as the helium flow rate increased from 6.8 ml/min 
to 20.3 ml/min at spaced ml/min intervals. With an increase 
of helium flow rate from 6.8 ml/min to 13.5 ml/min the 
increase in resolution was dependent on column length 
but not on the hydrocarbon number. For column lengths 
PeU mevcer to 5§0O meters the increase in resolution for 
n-nonane-decane, n-decane-hendecane and n-hendecane-do- 
decane was approximately 1.5 units. For column length 
O.5 meters the increase in resolution for the same series 
of compounds was approximately 1.0 unit. 

For an increase of helium flow rate from 13.5 ml/min 
to 20.3 ml/min the increase in resolution was dependent 
on column length and on the hydrocarbon number. For 
column lengths 1.0 meter to 5.0 meters the increase in 
resolution for n-nonane-decane, n-decane-hendecane was 


approximacely Ov7> units." The increase’ in>resolution 
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Fig. 4 Resolution vs programmed temperature rate when 
normalized net retention time for n-dodecane was fixed 
at 300 t 5 seconds. Columns from left to right are for 
eotunnog lengths tohtOes pel JO, 330, "and S7O meters “Rows 
from top, to-bottom are for n-nonane-decane, n-decane- 
hendecane and n-hendecane-dodecane. G, H and I are for 
helium flow rates of 20.3, 13.5 and 6.8 ml/min. 
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Fig. 5 Resolution vs column length when normalized net 
retention time for n-dodecane was fixed at 300 2 5 sec, 
Columns from left to right are for programmed temperature 
rates of 0.0 (isothermal), 2.0, 4.0, and’ 6.0 deg/min. 

Rows from top to bottom are for n-nonane-decane, n-decane- 
hendecane and n-hendecane-dodecane. M, N and O are for 
helium flow rates of 20.3, 13.5 and 6.8 ml/min. 
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Fig. 6 Same as Fig. 5 except columns from left to right 
are for programmed temperature rates of 8.0, 10.0, 12.0, 
and 13.0 deg/min. 
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for n-hendecane-dodecane was approximately 0.5 units. 

For column length 0.5 meters the increases in resolution 
for n-nonane-decane, n-decane-hendecane and n-hendecane- 
dodecane were approximately 0.4, 0.25 and 0.15 units. 

For a given pair of straight chain hydrocarbons and for 

a given increase in helium flow rate in the 6.8 ml/min 

to 20.3 ml/min range, the increase in resolution was inde- 
pendent of column length within the 1.0 to 5.0 meter range. 

For the separation of a given pair of straight chain 
hydrocarbons, for a given helium flow rate (6.8 - 20.3 
ml/min) and for a given programmed temperature rate (0.0 - 
13.0 deg/min) the resolution varied by less than one unit 
when the column length varied from approximately 1.5 meters 
to 5.0 meters. The resolution dropped sharply by several 
units when the column length decreased from approximately 
ie mecersrteo.10, 5 meters ive. swienr the *péencent liquid 
phase increased from approximately 19 to 50 percent. 

The optimum programmed temperature rate was dependent 
on which pair of straight chain hydrocarbons was being 
separated. As the straight chain hydrocarbon number in- 
creased the optimum programmed temperature rate decreased, 
For the separation of n-nonane-decane the resolution for 
all three helium flow rates increased gradually as the 
programmed temperature rate increased from 0.0 deg/min 


(isothermal) to 13.0 deg/min. By extrapolation the 
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optimum programmed temperature rate was approximately 

14 or 15 deg/min for column lengths of 0.5 meters to 3.0 
meters. The optimum rate was difficult to deduce for 

the 5.0 meter column. For the same range of’. programmed 
temperature rates (0.0 to 13.0 deg/min) and for the same 
pair of compounds (n-nonane-decane) the optimum column 
length shifted from 1.2 meters (isothermal) to 2.0 meters 
(13.0 deg/min) for the 20.3 ml/min helium flow rate and 
from 1.2 meters (isothermal) to 1.6 meters (13.0 deg/min) 
for the 13.5 ml/min helium flow rate. All of these column 
lengths are approximations. There was no shift in optimum 
column length for the 6.8 ml/min helium flow rate. The 
higher the helium flow rate the more was the shift in 
optimum column length. 

For the separation of n-decane-hendecane the reso- 
mrcion, Tomyalhethree: heddum <flowrratesand:doemradiiethe 
column lengths, increased gradually as the programmed 
temperature rate increased from 0.0 deg/min es ee aca 
to approximately 11.0 deg/min (optimum) and decreased 
slightly when the programmed temperature rate increased 
further from 11.0 deg/min to 13.0 deg/min. For the same 
range of programmed temperature rates (0.0 to 13.0 deg/ 
min) and for the same pair of compounds (n-decane-hen- 
decane), the optimum column length shifted from 1.6 meters 


(isothermal) to 2.5 meters (13.0 deg/min) for the 20.3 
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ml/min helium flow rate, 1.4 meters (isothermal) to 2.5 
meters (13.0 deg/min) for the 13.5 ml/min helium flow 
rate and 1.2 meters (isothermal) to 2.0 meters (13.0 deg/ 
min) for the 6.8 ml/min helium flow rate. All the above 
optimum column lengths are approximations. 

For the separation of n-hendecane-dodecane the 
resolution for all three helium flow rates in the column 
Jengthymange.0. 5emeters'ito 3.0emeters, wincreased gradually 
as the programmed temperature rate increased from 0.0 deg/ 
min (isothermal) to approximately 4.0 deg/min (optimum) 
and decreased when the programmed temperature rate in- 
ereased further from 4.0 deg/min to 13.0 deg/min. For 
the 5.0 meter column the optimum programmed temperature 
rate shifted from approximately 4.0 deg/min to 6.0 deg/min 
as the helium flow rate decreased from 20.3 ml/min to 6.8 
ml/min. For the same range of programmed temperature 
rates (0.0 to 13.0 deg/min) and for the same pair of 
compounds (n-hendecane-dodecane), the optimum column 
length shifted from 2.0 meters (isothermal) to 3.0 meters 
(13.0 deg/min) for the 20.3 ml/min helium flow rate, 1.8 
meters (isothermal) to 3.0 meters (13.0 deg/min) for the 
13.5 ml/min helium flow rate and 1.5 meters (isothermal) 
to beyond 5.0 meters (13.0 deg/min) for the 6.8 ml/min 
helium flow rate. All of the above optimum column lengths 


are also approximations. 
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4 Time Normalization With Net Retention Time For 
N-Dodecane Fixed at 800 + 10 Seconds 

Figs. 7 and 8 show the effect of helium flow rate, 
column length, isothermal temperature, and programmed 
remperature rare on the resolution of n=-nonane,t n-decane, 
n-hendecane, and n-dodecane when the net retention time 
for n-dodecane was fixed at 800 £ 10 seconds. 

The three helium flow rates used were again 6.8, 
13.5 and 20.3 ml/min. In all the cases shown the reso- 
lution=increased.as the helium flow rate increased from 
6.8 ml/min to 20.3 ml/min at spaced ml/min intervals. 

With an increase of helium flow rate from 6.8 ml/min 
to 13.5 ml/min the percent increase in resolution was 
agépendent on column length and ghydrocarbon number. For 
column lengths 1.0 meter to 5.0 meters the increases in 
resolution for _n=nonane-decaney n-decane=hendecane and 
n=hendecane-dodecane were approximately 2.0, 1.5 and 
O./5 units. For column lengths0.5 meters the increases 
in resolution for the same series of compounds were 
approximatvely* PRO sy "Os (5 "and O75 units. 

For an increase of helium flow rate from 13.5 ml/min 
to 20.3 ml/min the percent increase in resolution was 
also dependent on column length and hydrocarbon number. 
For column lengths 1.0 meter to 5,0 meters the increases 


in resolution for n-nonane-decane, n-decane-hendecane, 
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Fig. 7. Resolution vs programmed temperature rate when 
normalized net retention time for n-dodecane was fixed 
at 800 = 10 seconds. Columns from left to right are for 
Covunntbensuns -ofr70T>, TsO, 5.0, and 570 meters. Rows 
from top to bottom are for n-nonane-decane, n-decane- 
hen@ecane “and *n-hnendecane-—dodecane. ~ J, K°and: Lvare for 
Hetitm/fiow rates of 20.3, 13.5 and 6.8 ml/min. 
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Fig. 8 Resolution vs column length when normalized net 
retention time for n-dodecane was fixed at 800 = 10 sec. 
Columns from left to right are for 0.0 (isothermal), 2.0, 
HO, .and 6.0 deg/min. Rows from top to bottom are for 
n-nonane-decane, n-decane-hendecane and n-hendecane-dodecane. 
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and n-hendecane-dodecane were approximately 1.0, O.4 
and 0.25 units. For column length 0.5 meters the in- 
crease in resolution for n-nonane-decane was approxi- 
mately 0.2 units and was insignificant for n-decane- 
hendecane and n-hendecane-dodecane. For a given in- 
crease in helium flow rate in the 6.8 ml/min to 20.3 ml1/ 
min range, the increase in resolution became more dramatic 
as the straight chain hydrocarbon number decreased. For 
argiven pair of ‘straight chain hydrocarbons.,and ifor a 
given increase in helium flow rate in the 6.8 ml/min to 
20.3 ml/min range, the percent increase in resolution | 
was independent of column length within the 1.0 meter to 
5.0 meter range. 

For the separation of a given pair of straight 
chain hydrocarbons, for a given helium flow rate (6.8 to 
20.3 ml/min) and for a given programmed temperature 
rate (0.0 to 6.0 deg/min) the resolution varied by less 
than one unit when the column length varied from approxi- 
mately 1.5 meters to 5.0 meters, (The resolution dropped 
sharply by several units when the column length decreased 
from’ approximately 175° meters to Orcd5ameters: ijveca whens the 
percent liquid phase increased from approximately 19 to 
50 percent. 

The optimum programmed temperature rate was dependent 


on which pair of straight chain hydrocarbons was being 
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separated. As the straight chain hydrocarbon number in- 
creased the optimum programmed temperature rate decreased. 
For the separation of n-nonane-decane the resolution 

for all three helium flow rates and all the column lengths 
increased with programmed temperature rate from 0.0 deg/ 
min (isothermal) to 6.0 deg/min (optimum). For the same 
range of programmed temperature rates (0.0 to 6.0 deg/min) 
and for the same pair of compounds (n-nonane-decane) the 
optimum column length shifted from 2.3 meters (isothermal) 
to beyond 5.0 meters (6.0 deg/min) for the 20.3 ml/min 
helium flow rate, 2.0 meters (isothermal) to beyond 5.0 
meters (6.0 deg/min) for the 13.5 ml/min helium flow rate 
and 1.5 meters (isothermal) to 4.0 meters (6.0 deg/min) 
for the 6.8 ml/min helium flow rate. All of these column 
lengths are approximationsy: Closecobservation of-the 
graphs showed that the higher the helium flow rate the 
more was the shift in optimum column length. 

For the separation of n-decane-hendecane the reso- 
Luvlon Cor ‘alle the: hektiumetilow rates and 3all the column 
lengths, increased as the programmed temperature rate 
increased from 0.0 deg/min (isothermal) to 4.0 deg/min 
(Optimum) and decreased when the programmed temperature 
rate increased further from 4.0 deg/min to 6.0 deg/min. 
For the same range of programmed temperature rates (0.0 


to 6.0 deg/min) and for the same pair of compounds 
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(n-decane-hendecane), the optimum column length shifted 
from 2.8 meters (isothermal) to 5.0 meters (6.0 deg/min) 
for the 20.3 ml/min helium flow rate, 2.8 meters (iso- 
thermal) to 4.0 meters (6.0 deg/min) for the 13.5 ml/min 
helium flow rate and 2.5 meters (isothermal) to 3.5 meters 
(6.0 deg/min) for the 6.8 ml/min helium flow rate. All 
the above optimum column lengths are approximations. 

For the separation of n-hendecane-dodecane the reso- 
Jution,for,all, three hediumsaf Loweratesrandeallathe, column 
lengths, increased as the programmed temperature rate in- 
creased from 0.0 deg/min (isothermal) to 2.0 deg/min (op- 
timum) and decreased when the programmed temperature rate 
increased further from 2.0 deg/min to 6.0 deg/min. For 
the same range of programmed temperature rates (0.0 to 
6.0 deg/min)and for the same pair of compounds (n-hende- 
cane-dodecane), the optimum column length remained fixed 
at approximately 3.0 meters for the 20.3 and 13.5 ml/min 
helium flow rate and at approximately 4.0 meters for the 


6.8 ml/min helium flow rate. 
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5 Time Normalization With Net Retention Time For 
N-Dodecane Fixed at 1600 * 15 Seconds 

Figs. 9 and 10 show the effect of helium flow rate, 
column length, isothermal temperature, and programmed 
temperature rate on the resolution of n-nonane, n-decane, 
n-hendecane, and n-dodecane when the net retention time 
for n-dodecane was fixed at 1600 = 15 seconds. Fig. 9 
is. incomplete because as pointed out in the introduction 
initial column temperatures were restricted to those 
above room temperature. 

The three helium flow rates used were again 6.8, 
13.5 and 20.3 ml/min. In all the cases shown the reso- 
lution increased with helium flow rate from 6.8 ml/min 
to 20.3 ml/min at spaced ml/min intervals. With an increase 
in helium flow rate from 6.8 ml/min to 13.5 ml/min the 
percent increase in resolution was dependent on.column 
length and the hydrocarbon number. For column lengths 
1.0 meter to 5.0 meters the increases in resolution for 
n-nonane-decane, n-decane-hendecane, and n-hendecane- 
dodecane, were approximately 2.0, 1.25 and 0.5 units. 
For column Jeneth 0.5 meters the ancreases in resolution 
for the same series of compounds were approximately 1.0, 
O27 5-ana"0.5"units: 

For an increase of helium flow rate from 13.5 ml/min 
to 20.3 ml/min the percent increase in resolution was 


also dependent on column length and the hydrocarbon number. 
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Fig. 9 Resolution vs programmed temperature rate when 
normalized net retention time for n-dodecane was fixed 

at 1600 + 15 seconds. Columns from left to right are 

Por cosimn Leneuis OL” Ovo, sO > 320) rand, UNMeters. 

Rows from top to bottom are for n-nonane-decane, n-decane- 
hendecane and n-hendecane-dodecane. S, T and U are for 
helium flow rates of 20.3, 13.5 and 6.8 ml/min. 
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0.0 DEG/MIN 2.0 


Resolution 


1 3 5 1 3 5 
Column Length (meters) 


PaieoxriOatkesolutiontivescolumelengthawhen normalized net 
retention time for n-dodecane was fixed at 1600 = 15 sec. 
Columns from left to right are for 0.0 (isothermal) and 
2.0 deg/min. Rows from top to bottom are for n-nonane- 
decane, n-decane-hendecane and n-hendecane-dodecane. 

V¥, Wand X are for helium. fiow rates of 20.3, 13.5. and 
6.8 ml/min. 
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For column lengths 1.0 meter to 5.0 meters, the increases 
in resolution for n-nonane-decane, n-decane-hendecane 

and n-hendecane-dodecane were approximately "0.75;) O35 tand 
0.2 units. For column length O.5 meters the increase 

in resolution for the same series of compounds was insig- 
nificant. For a given increase in helium flow rate in 

the 6.8 ml/min to 20.3 ml/min range, the increase in reso- 
lution became more dramatic as the hydrocarbon number 
decreased, For a given pair of straight chain hydrocarbons 
and for a given increase in helium flow rate in the GS 
ml/min to 20.3 ml/min range, the percent increase in 
resolution was independent of column length within the 

1.0 meter to 5.0 meter range, 

For the separation of a given pair of straight chain 
hydrocarbons, for a given helium flow rate (6.8 to 20.3 
ml/min) and for a given programmed temperature rate (0.0 to 
2.0 deg/min), the resolution varied by less than one unit 
when the column length varied from approximately 1.5 
meters to 5.0 meters. The resolution dropped sharply 
by several units when the colum length decreased from 
approximately 1.5 meters to 0.5 meters i.e. when the 
percent liquid phase increased from approximately 19 to 
50 percent. 

For the separation of n-nonane-decane the resolution 


for all three helium flow rates and all the column lengths 
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increased as the programmed temperature rate increased 
from 0.0 deg/min (isothermal) to 2.0 deg/min. By extra- 
polation the optimum programmed rate for the 6.8 ml/min 
helium flow rate (column length 3.0 meters to 5.0 meters) 
was approximately 4.0 deg/min. The optimum programmed 
rate for the other two helium flow rates could not be 
determined, For the same range of programmed temperature 
rates (0.0 to 2.0 deg/min) and for the same pair of com- 
pounds (n-nonane-decane) the optimum column length shifted 
from 3.5 meters (isothermal) to beyond 5.0 meters (2.0 
deg/min) for the 20.3 ml/min helium flow rate, 2.5 meters 
(isothermal) to 5.0 meters (2.0 deg/min) for the 13.5 
ml/min helium flow rate and 2.0 meters (isothermal) to 
5.0 meters (2.0 deg/min) for the 6.8 ml/min helium flow 
rate. All of these column lengths are approximations. 
Here again the higher the helium flow rate the more the 
shift in optimum column length. 

For the separation of n-decane-hendecane the reso- 
tution Tor allthe helium flow rates and all tthe’ column 
lengths increased as the programmed temperature rate 
increased from 0.0 deg/min (isothermal ) to 2.0 deg/min. 
The optimum programmed temperature rate for the 6.8 ml1/ 
min helium flow rate (column length 3.0 meters to 5.0 
meters) was approximately 2.0 deg/min. The optimum 
programmed temperature rate for the other two helium 


flow rates could not be determined. For the same range 
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of programmed temperature rates (0.0 to 2.0 deg/min) and 
for the same pair of compounds (n-decane-hendecane) the 
optimum column length shifted from 4.0 meters (isothermal) 
to 5.0 meters (2.0 deg/min) for the 20.3 ml/min helium 
flow rate, there was. no significant, shift in optimum 
column length for the 13.5 ml/min helium flow rate and 
shifted from 3.5 meters (isothermal) to 4.0 meters (2.0 
deg/min) for the 6.8 ml/min helium flow rate. All the 
above optimum column lengths are again approximations. 
For the separation of n-hendecane-dodecane the reso- 
bucion.tor,ai) .boaeece helium flow rates, and .all the .column 
lengths, increased as the programmed temperature rate 
increased from 0.0 deg/min (isothermal) to 2.0 deg/min. 
The optimum programmed temperature rate for the 6.8 ml/ 
min helium flow rate (column length 3.0 meters to 5.0 
meters) was approximately 2.0 deg/min. The optimum pro- 
grammed temperature rate for the other two helium flow 
rates could not be determined. For the same range of 
programmed temperature rates (0.0 to 2.0 deg/min) and 
for the same pair of compounds, there was no shift in the 
optimum column length when the programmed temperature 
rate increased from 0.0 deg/min (isothermal) to 2.0 deg/ 


min for all three helium flow rates, 
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6 Comparison of Experimental and Theoretical Curves 
Figure 11 shows the experimental and theoretical 
resolution between n-hendecane and n-dodecane versus 
column length for three (300, 800 and 1600 seconds) nor- 
malized net retention times and three (6.8, 13.5 and 20.3 
ml/min) helium flow rates for isothermal column tempera- 


tures<.,.R.is...the..resolution..calculated-from the equation 


R = (tpo-tp,)/0.5(w,+w,) Egy uit 
where ty and W, are the net retention time and peak width 
for the first peak and Cro and W, are the net retention 


time and peak width for the second peak (5). Rl is the 
resolution calculated from the equation 

1 #17) {1/H, (a-1/a) (kp/K5+1) Eq. 2 
the parameters of which are defined in the introduction. 


Rew is the theoretical resolution calculated from the 


equation derived by Grushka, Yepes-Baraya and Cooke (3) 
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another theoretical resolution Ryo where Ryo = 
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Fig. 11 Resolution between n-hendecane and n-dodecane 

for normalized net retention times of 300, 800 and 1600 
seconds and helium flow rates of 6.8, 13.5 and 20.3 ml/min. 
Poe equapions for the tour curvessare Roe Bquation.'h, 


Rl = Equation 2, Ray, = Equation 3, and Ray. = Equation Ne 


All the curves are for isothermal conditions. 
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which gives a curve more related to the experimental 
results of this work, This equation was obtained by 
tenet and error using the second computer program in the 
appendix. 


The curve Rea im Figure (lls ai appreximacve <strareny 


line the slope of which increases with normalized net 
retention time... The experimental resolution R or Ri 

on the other hand increases sharply with column length 

to about 1.5 meters and then levels off and decreases 
Slightly with increasing column length. It was decided 

to add a third term to the denominator in Equation 3 which 
would increase the resolution at short column length and 
decrease resolution at longer column length. This was 
first partly overcome by making the second term in the 


denominator of Equation 3 less dependent on Kn since Ky 


decreases sharply with increasing column length (see Tables 


1 and 2). The factors k,(k,+1) in the second term were re- 


2 3 


placed by foe Since k, changes so rapidly with co- 


2 


lumi Length, 10 was decided that Ky Should be an important 


factor in the third term. To minimize the rapid decrease 


of Ky with increasing column length the third term took 


the form of a ratio k/K, each k, to a different power. 


2 
Since the experimental resolution R or Rl is also dependent 


on the normalized net retention time it was decided that 


the Ky in the numerator of the third term should be to 
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some power involving Cpo- To minimize the large numerical 
change in tro» this power also took the form of a ratio 
tpo/tpo each tp, to a different power if necessary. The 
final form of this ratio was (eee ee Vitae. To intro- 
duce :a curvature into the curve in the vicinity of écodumn 
length of 1.5 meters it was decided to use: a log factor 
which was a, function ol column, jlensth.~ The cfinal forni 
of this factor was log(1+L/15). Since the experimental 
PesoluciLon is alesocaliuncuions or heliumitlow rate, this 
was best overcome by multiplying the first term in the 
denominator of Equation 3 by the compressibility cor- 
rection factor j. Thus the end result was Equation 4. 
Tie (ODN LOUS Chau there leo. ouL i. a Ou, Of “POOm Lor wine 
provement. Referring to Figure ll, Rayo increases too 
rapidly with increasing helium flow rate at long norma- 
zed net retention times. This leads one. to speculate 
that the third term (or even the first or second term 
or possibly a fourth term) in the denominator of Equation 
4 should contain yet another factor e.g. the compressi- 
bidity. factor, jnasnactunction of normalized net uretention 
time toe 

The relevant data used to calculate the curves in 
Fig. 11 are shown in Tables 1 and 2. The constants which 
were also used are S = 14500. om-/em, Ps = 0.240 gm/cc, 
M, = MB My = Tn ot, Eee A Too = 659.0" °K, 


Vol = 57.8 cc/g mole, and Voo = 718. cc/g mole. 
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The relative retention ratio a was calculated from 
the net retention times of n-dodecane and n-hendecane 


ee oot The compressibility factor j was 


Ro/*R1: 
calculated from the equation j = 3/2(((P,/P,)°-1)/ 
((P,/P,)>-1)) where P, is the outlet column pressure and 


P, is the inlet column pressure (6). The inlet column 


al 
pressure was measured with a Matheson P/N 63-4111 pressure 
gauge to which had been attached a pressure line and 
needle. This was inserted through the rubber septum. 


The diffusion coefficient of a solute in the gas 


phase was calculated from the formula (7) 


0.43(2/100)**92( (1/1, )+(1/M) )°* 
0.1405 ((y 


)O-4 


O.4,2 
P(T 4 To5/10000 ) 1/100 + Vgo/100) ) 


wnere T is the absolute temperature in °K, P is the 


column outlet pressure (8), M, and M, are the molecular 


weights, Tol and Too are, the critical temperatures in: °K, 


and V,, and V,, are the critical volumes in em?/g mole 


al 
of: the eluent gas and diffusing solute. 


The value of DS was obtained by using the assumption 


at 
value of 0.6, 


that D, = 10°°D, The obstruction factor y was given a 


When the dead volumes i.e. the retention volume of 
the air peak for the three helium flow rates were extra- 
polated to zero column length the new dead volumes were 
6.01, 5.62 and 4.90 mls (rétention times 53.0, 25.0 and 


14.5 seconds) for the 6.8, 13.5 and 20.3 ml/min helium 
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flow rates. These new dead volumes represented the total 
dead volume of the gas chromatographic system minus the 
dead volume of the columns. The formulas used to calculate 
the capacity ratio k, for the GYo siise 5 sand. Os samy iar? 
helium flow rates were 

SUN Ss SONA emis aeN k 

WEN 25. Oy / (Va -25.0)) and 

gee Vg 14.5) / (V1 5) 


where V is the retention time (seconds) of n-dodecane and 


a 
| 


Ws 
Il 


V,. is the retention time (seconds) of the air peak. 
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( General Conclusions When Using Programmed Temperature 
Time Normalization 


1) Resolution increased with helium flow rate. This 
observation was independent of the normalized retention 
time, column length and hydrocarbon number. 

2) For a given pair of straight chain hydrocarbons and 
for a given increase in helium flow rate in the 6.8 to 
OD rew ml/min range,.the, percent increase in resolution 
was the. same.for.column lengths.from.1.0.to.5.0) meters 
DUT WasmOlllernent fom the, 0.5. merer, column .lengrh nas 
observation was also independent of the normalized re- 
Lenton, ime. | 

3) For an increase in helium flow rate from 6.8 to 
13.5 ml/min the percent increase in resolution was 
dependent on column length from 0.5 to 1.0 meter. The 
increase in resolution was also dependent on the hydro- 
carbon number for net retention times of 800 + 10 and 
1600 + 15 seconds but not for 300 + 5 seconds. 

4) For an increase in helium flow rate from 13.5 to 
20.3 ml/min the percent increase in resolution was 
dependent on column length from 0.5 to 1.0 meter and on 
the. yhydrocarbon.number. This was the case for all net 
retention times from 300 + 5 seconds to 1600 * 15 seconds. 
5) For a given increase in helium flow rate in the 6.8 
GO..20.. 3 ml/min range, the percent increase in resolution 


became more dramatic as the hydrocarbon number decreased, 
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6) For the separation of a given pair of straight chain 
hydrocarbons, for a given helium flow rate (6.8 to 20.3 
ml/min) and for a given programmed temperature rate (0.0 
to 13.0 deg/min, 0.0 to 6.0 deg/min and 0.0 to 2.0 deg/min 
for normalized net retention times of 300 ¢ 5, 800 1 
and 1600 + 15 seconds), the resolution varied by less 

than one unit when the column length varied from Seat 
mately 1.5 to 5.0 meters. 

7) For normalized net retention times of 300 + 5 to 

1600 + 15 seconds, the resolution dropped sharply by 
several units when the column length decreased from 
approximately 1.5 to 0.5 meters i.e. when the percent 
liquid-phase increased from approximately: 19 to 50-percent. 
8) The optimum programmed temperature rate was dependent 
on which pair of compounds was being examined. As the 
straight chain hydrocarbon number increased the optimum 
programmed temperature rate decreased. Also, for a ‘given 
pair of hydrocarbons, the optimum programmed temperature 
rate decreased as the net retention time increased up to 
about 1600 seconds; here there was indication that the 
limiting optimum programmed temperature rate was approxi- 
mately 2.0 deg/min. 

9) Generally the higher the helium flow rate the more 
was the shift to higher values in the optimum column 


length. 
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Computer program for the calculation of HETP and resolution 
Prologue 

The computer program written in Fortran can be divided 
into three major Do loops. The first loop incorporates 
the entire programm except for the first few cards. The 
second loop is within the first loop and the third loop 
is within the second. 

| The data for this program was obtained from a printer 
which was connected to a digital integrator which in turn 
was connected to the detector through a Gow-Mac power 
supply. A change of 1 mv in the detector would give a 
response of 1000 counts on the integrator. The system 
was set up so that the integrator would integrate for 
aj weconds print ioubrchetotalinumber tofiicountspiweset the 
integrator to zero and repeat the process. In this way 
a numerical outline of a chromatogram could be obtained. 
If the timer (built into the integrator) was started at 
the point of injection of a sample, the time (in seconds) 
as well as the number of counts could both be utilized 
by the computer program. 

Two cards are read into the computer before the first 
DOALOOp.. 4-Card “bicon tains “tie ol Lowine int ormatons 
Columns 1-5 contain an integer which indicates how many 
data numbers will appear on each card when the base line 
and numerical outline of the chromatogram are read into 


the computer. For this program this number can be either 
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8,15,0r 165 The ‘base ‘line rand» peak outline will be 
mentioned again later. 

Columns 6-10 contain an integer which indicates the number 
of sets of data (one set of data consists of three injections) 
that will be read into the computer. This will govern the 
number of times the first major Do loop will be repeated. 
The next four numbers must be real numbers. 

Columns 11-20 contain the helium flow rate (ml/min). 
Columns 21-30 contain the detector current (ma). 

Columns 31-40 contain the injection port temperature (°C). 
Columns 41-50 contain the sample size (m1). 

Columns 1-48 on the second card to be read in will 
contain the names of the sample components (a maximum of 
of 4 components in one computer run) that are being se- 
parated. The first component is written in columns 1-12, 
the second in columns 13-24, the third in columns 25-36, 
and the fourth in columns 37-48. 

Inside the first major Do loop one card is read: into 
the computer. This card contains the following information. 
Columns 1-5, 6-10 and 11-15 contain integers of the date 
(day, month and year). 

Columns 16-20 contain either the integer 1 or the inte- 
ger 2. If it is 1 the the data was obtained under iso- 
thermal conditions. If it is 2 then the data was obtained 
under programmed temperature conditions and the isothermal 


temperature (columns 26-35) read into the computer is 
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taken as the initial temperature of the column. 

Podumnbsie lees contain el thenmets the> integers! 2388 ox ie 

If it is 2 then the data for only two peaks (two components) 
Wadinibe read intortheicomputer ..edf. itCisagithen) the data 
for three peaks (three components) will be read into the 
peapureta If 4 then the data for four peaks (four com- 
ponents) will be read into the computer. The next five 
numbers must be real numbers. 

eo inimnel os S5 contain the isothermal temperature (°C) of 
the coun: 

Columns 36-45 contain the programmed temperature rate (deg/ 
min). 

Columns 46-55 contain the sensitivity of the detector. 
Columns 56-65 contain the column length (m) and 

Columns 66-75 contain the percent liquid loading in the 
column. 

Inside the second major Do loop three types of data 
are read into the computer. Columns 1-5 on the first card 
contain a real number indicating the time (in seconds) of 
the air peak. Columns 1-5 on the second card contain an 
integer indicating the number of baseline values (counts) 
that will be read into the computer. The next card (one 
or more depending on how much baseline is read into the 
computer) contains the baseline values. The number of 


values per card has been previously determined. 
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Inside the third Do loop four types, of data are 
read into the computer. Columns 1-5 on the first card 
contain an integer indicating the time in seconds at the 
beginning of a chromatographic peak. Columns 1-5 on the 
second card contain an integer indicating the number of 
peak values that will be read into the computer. One 
peak value corresponds to one print out of the printer. 
Columns 1-5 on the third card contain a real number the 
value of which depends on a combination of several factors. 
The final value of this number depends on the print rate 
of the printer and on whether every single one or every 
SeconarOPeuiir Ol 1 Our ul CLC, OL Ulle. prin OuUL Values 
of the peak are read into the computer, For example if 
the print rate ofthe printer was 1 print per 2 seconds 
and every single value of the peak print out was read 
into the computer, the”value of ‘this “number~wotuldzbe 2, 
i.e. the total number of print out values from the be- 
ginning of the peak to the peak maximum multiplied by 2 
and added to the time at the beginning of the peak would 
give the retention time of that peak. If the peak was 
very broad so that only every tenth peak value (with 
print rate = 1 print/2 seconds) need be read into the 
computer ‘then the value of this “number *would “be 2O. The 
fourth card (one or more) contains the peak values of the 
chromatogram. The number of peak values per card must be 


the same as for the baseline values, 
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The number of times the third Do loop repeats itself 
depends on the number of peaks or components (maximum of 4) 
for each injection. If 4 peaks are to be read into the 
computer then the third Do loop will repeat 4 times. The 
second Do loop will repeat itself 3 times (for 3 injections). 
The number of times the first Do loop will repeat depends 
on the number of sets (each of 3 injections) that will 


be read into the computer. 


Epilogue 


The program prints out all the data read into the 
Computer... in, additions it will calculate’ and print’ out 
for each injection and also the average of the three in- 
jections the following: relative peak areas, retention 
time of each peak, peak heights, baseline peak width, 
resolution using the formula R = (tpo-tp,)/0.5(wi+W,) 
theoretical plates, and peak skew. It prints out the 
average retention time of the air peak, average net re- 
tention time of each peak, average capacity ratio using 
the formula k = (V-V55)Vas: HETP using the formula 
HETP = 16(t,/w)*, and effective plate number. In ad- 
dition it prints out miscellaneous information such as 
the positive and negative slopes of the peaks and the 
(x,y) coordinates on the peaks used to obtain the slopes. 


This, information proves userul in tracking down the source 


of any errors that may appear in the output. 


ae 


ae 


abe 


tloedt atesqet qoc i om eo 


to mumixam) a JroH0qHo® * 


Soot ed og on 


* 1 


tw qool of sant? ist 


= 
a 


& ai a - "} 
(amottoetat & Yo fiuat 3 ; : 


steb oat ifs tuo ad cikag MBTR 2 i 3 
Hsiuolso Liiw HS months at. 
sdt+ oate bee nrottost mt a 
tsiex :golwolfoet edd ade 
tkiepeda .&sra lsc Raq ise q oes a 
[ eer : « cer i 
t)} = A siumitel sod en tew nic obit ule 
JI. wenle neq bne ,petelig Deotde 
M. 


yr sk a as a o 5 ety * ry iT att nokineten og - 


s2eatevs ,teeq dose to saitd mob 


ola! evpsasiieorkm duo hair 
eit lo saqels evideaen bas, ov id 


ot beew eussq sit no aadsaths 


gitdostt ot Lyuteay eeveg moldemuot 


“id ot ta9qgqs Yam dade a 


ei 


REAL X (7000) ,AREA(12) ,AVGE(4), BASE(48) , RET (12) ,PEAK(12) ,M(24) 

C,WWOCT (3) ,WWNON (3) ,WWDEC (3) , WWHEN(3) ,RRON(3) ,RRND(3) ,RRDH (3), 
CTEPOCTY3) ,THFNON (3) , THPDEC (3) , THPHEN (3) , TAIR (3), AM(8) « AW (8), 
CTA ,SKOCT (3) ,SKNON (3) , SKDEC(3) , SKHEN (3) ,AA (16) ,BB (8) ,XXX{4,12), 
CYYY(2,12),TRE(4) ,KPFA (4) ,WONDH (4) ,ARONDH (4) ,CPONDH(4),SKONDH({4), 
CHECNDH (4) ,EPONDH (4), TPONDH (4) ,W(24) 

INTEGER T,CARD (20) 

READ (5, 20 Cy NEM NOS g HER DC, rl, oS 

200 FORMAT (215, 4F10. 2) 
Hero, tC tT) s(CARD CL itu 15.12) 
aoe 4 FORMATO RS 
DO*900 “TA=1,NS 
READ (S520 1) COA 1, TOD Asy LDAS, KR KKE CT, PRA,SE,CL,PLP 
201 TORRAT (SES oF 10.2) 
GO T0 (985,986) ,KK 
$86 WRITE (6,987) TDA, IDA2,1DA3,PRA,CT 
GA) FORMAL toe DATE=", 2X,12, 1h, i2, ike lope ks MPROGRAMNMED. TEMP RATE=", 
C2k 5 F5. ee DEGREES CENTEGRADE PER MINUTE WITH INITIAL TEMPEH", F7.2 
Cycks  VEoR Peo CUNT IGRADE ©) 
WRITE (6,988)CL,SE 
peo FORMAT OC OuUNN PENG TH="Foe2 tap EET ERS (92k ,{OSENSLTIVITY =", 75. 2) 
GO TO 992 
265 WRITE (fo, USS IDAN, UDAZ, PDAS CT, CL, SE 
Goo FORMAT Ws F hes 2k, le, bk log (kpc eke LOOTHERMAL COLOMN TEMP=', 
Chi. 1, thy LEGRGEES CENTIGRADEN 2 ¥, ‘COLUMN LENGTH=",F/.2, 1X, "METERS", 
C2, Vo ENS TEIN Ty = 995.2) 
32) Wr Ee fo 990) aie ne DC PS 
99) FORA TP MeL PLOW RAD He vole ge Xs (MLS PER MINUTE*?, 2%, * DETECTOR 
CCURRENT=EC SPO. 1a 2k, UM LLLIAMPS!, 2X, "ENJECTION PORT TEMNP=*, Fb. 1, 2X, 
CULDPGRENS CENTAGRADE' ) 

WRITE(6,991) SS, (CARD (I) ,I=1, 12) 

991. FORMAT ("SAMPLE SIZE INJECTED WAS*,F6.2,2X;,'MICROLITER?, 2X, 12A4) 

WRITE(6,912) PLP 

S12 ORMAT (8 Tt of COLUMN Toe rue sy 2k Ge CENT UCH9S ON. CHROMOSORS 
CWAW 60-80 MESH IN 4VYS" TNCH’CU TUBING") 
915 NI=0 

NH=0 

IP=0 

Is=0 

DO 800 IN=1-,,3 

RERD({5, 202) TA 

Z02 FORMAT FS. 1) 

READ (5,203) TTA 

203 FORMAT (15) 
LF (NFS EO. 169° GO TO 1502 
TF(NFM. EQ. 15)* GO TO'1503 
READ(5, 1504) (BASE(I) ,I=1,ITA) 
1504 FCEMAT (8 F10.0) 

GOS 70 4505 
1502 READ{S. 1506) (BASE (1) 52-1, FTA) 
1506 FOEMAT(16F5.0) 

GOLeO FSC" 

1503 READ(5, 1507) (BASE(I) ,1=1,ITA) 
1507 FPORUATAC! SPS. 0) 
1505 TRIR(I8y TA 
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SBASE=0.0 

IT=IN 

DO 801 I=1,ITA 

SBASE=SBASE+BASE (TI) 

AVEBASE=SBASE/ITA 

WRITE(6,802) IN 

FOTMAT (‘-',*THE BASELINE VALUES (CCUNTS) FOR INJECTION',I3) 
WRITE (6,803) (BASE(I) ,I=1,ITA) 

FORMAT (1X,2CF6. 0) 

WPITE (6,804) AVBASE 

FORMAT('THE BASELINE AVERAGE =!,F7.0,2X,"SECONDS') 
WRITE (6,805) TA 

FORMAT (‘THE AIR PEAK RETENTION TIME =',F7.2,2X,'SECONDS') 
DO 700 IM=1,KKK 

READ(5,72C)T 

FORMAT (15) 

READ{5,721)N 

FORMAT (15) 

RFAD.(5, 722) DF 

FORMAT (F5. 2) 

NNN=T4N-1 

IF (NFM.EQ.16) GO TO 1510 

IF (NFM.EQ.15) GO TO 1511 

READ (5,1512) (X(I) , I=T, NNN) 

FORMAT (8F10.0) 

GO TO 1513 

READ (5, 1520) (X(I) ,I=T,NNN) 

FORMAT (16F5.0) 

GO TO 1513 

READ(5,1514) (X(1) ,I=T, NNN) 

FORMAT (15F5.C) 

Go m0l210,211, 212,213) ,18 
WRITE(6,214) (CARD(I) ,I=1,3) ,IT 

FORMAT ('04,72,'THE PEAK VALUES (COUNTS) FOR',T32,3A4,T45, 


CUrnuvecTr ion! , 13) 


Gol TH! 222 
WRITE (6,215) (CARD (I) ,I=4,6),1T 
FORMAT ('O', 12,'THE PEAK VALUES (COUNTS) FOR',1T32,3A4,T45, 


C'INJECTION*,13) 


Go MO) 209 
WRITE (6,216) (CARD(I) ,I=7,9),1T 
FORMAT ('0',172,'THE PEAK VALUES (COUNTS) FOR',T32,3A4,T45, 


C*INJECTION',13) 


GO! To’ 2212 
WRITE(6,218) (CARD(I) ,I=10,12) ,IT 
FORMAT (*0',12,*THE PEAK VALUES (COUNTS) FOR?,1T32,3A4,T45, 


C'INJECTION',13) 


WRITE(6,219) (X(I) ,I=T,NNN) 

FORMAT (1X, 17F7.0) 

WRITE(6,220)N 

FORMAT ("NUMBER OF CURVE DATA POINTS READ IN=',1I5) 
WRITE(6,221)T 

FORMAT ('TIME AT BEGINNING OF PEAK=",1I5,2X,'SECONDS') 
KE=AVEASE 

K B=0 
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DO 704 T=T,NNN 

DER oe XC) ) Ot et OS 
KB=KB+1 

DEH Bee Ves) eGCl TO 705 

GO TO 704 

KR=0 

GO TO 704 

KB=0 

TPUxAI tte GiakKE) GO TG 706 
GO TO 704 

IP=I+1 

KE=X (I +1) 

CONTINUE 

NI=WI+1 

PEAK (NI) =X (IP) -AVBASE 
AR=C.0 

NN=T+N-1 

Sees GOs Ube bos Je 
Lela. Oe 0) GU! LO vero 

DO 707 I=T,NNN 

AR=AR+X (I) 

DO 710 I=T,NN 

ARS ARS CCX AT) +X (14 1) ) 7 2..0) FD E-160) 
AR=AR-AVEASF*N*DF 

GO TO. 711 

DO 714 I=T,NNN 

AR=AR+X (1) 

AR=AR-AVEASE*¥N 

GO TCO 711 

NBE=0 

DOL iS LT NNN, 2 

NE=NB+1 

AKR=AR+X (T) 

AR=AR-~AVEBASE*¥NB 

NH=NH+1 

AREA (NH)=AR 

PK 1=PEAK (NI) *0.507 
PK2=PEAK (NI) *0.707 

IR=IR+1 

DG. 24 °° l=T1P 
IF{PK1.GT.({X{2I) -AVBASE)) GG TO 724 
PRISIPRI {x (1-1) AV BAS E)) ad 021) -AVSAS ©) — (% 01-1) —AVBASE)) 
YYY (1,18) =PK1 

Rael s (l=) TeP ne De 

KXX (1,18) =xXxX1 

GU TO, 519 

CONTINUE 

Von de l= ty Le 

LEI PEKZ. Gb)» (Xx (LY—AV BASS) ))° GO TO 725 
PH2= (PK 2-({X (I-1) -AVBASE) ) 7 ((X (I) -AVBASE) -—(X (1-1) -AVBASE) ) 
YYY(2,1R)=EK2 
XX2=T+({ (1-1) -T+PH2) *DF 

XK i 2, LB) HXKZ 

GO. TO 544 


725 CONTINUE 
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DO 726 I=IP,NNN 

AEYPK2.6Ts.(X (1)=—AVBASE)) GO TO 727 
G@ TO 726 

PHS= (PK2— (X (I) ~AVBASE) ) / ( (X (1-1) ~AVBASE) — {X({1) -AVBASBE) ) 
RXS=T+ ((1-—1) -L+1. O- PHS) #DF 

XXX (3,1R)=XX3 

GO"TO 512 

CONTINUE 

DC 728 I=IP,NNN 

OF (PRISGTS UXT) -AVBASE)) <GO TO 729 
GO TO 728 

PH4= (PK1- (X (I) ~-AVBASE)) /((X (1-1) ~AVBASE) - (X (1) ~AVBASE) ) 
XX4=T+ ( (1-1) -T+1.0-PH4U) *DF 

XXX (4,18) =XxX4 

GO" POPS 3 

CONTINUE 

IS=1S+1 

M(TS)=(PK2-PK1)/(XX2-XX1) 

WAITS) = OC LS) KX I-PK 1) JB (LS) 
TS=15+1 

M (IS) = (PK 1-PK2) / (XX4-XX 3) 

WCTS) = (4% (TS) *XX3-PK2) /K (IS) 

Peo NG) = Cs S 1) FX KX 21S) KX 3)7 (8 Gs t-te (1S)) 
CONTINUE 

CONTINUE 

A= (LAIR (1) #TAIR (2) +TAITR{(3)) 73.9 

CO WO (223,,223,224,,225) ,, KKK 

JB=5 

JE=9 

SOeTO W225 

JB=4 

JE=7 

GO" T0-1225 

JB=3 

JE=5 

NE=0 

DO 1224 I=1,KKK 

TRE (LI) = (RET (NK*+1) t+RET (NK+JB) +RET (NKtJE) ) /3.0 
NK=NK+1 

Gin tor (2:27, 227,228,229) KKK 

JB=5 

JF=9 

Go eu M229 

JB=4 

JF=7 

GO 7TO-. 1229 

JB=3 

JE=5 

NK=0 

DO 1228 I=1,KKK 

KPEA (1) = (PEAK (NK+1) *PEAK (NK+JB) + PEAK (NK+JE) ) /3.0 
NK=NK+1 

NA=0 

NK=0 

GO TO (231, 231,232,233) , KKK 
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JB=5 

JE=9 

GOs TO, 524 

JR=4 

JF=7 

GO TO 524 

JB=3 

JE=5 

DO, S202 Tedy ih 

DO 521 J=1,KKK 

NA=NA+1 

AA (NA) = (XXX (I, NK4+1) +XXX (L,NKt3B) +XXX (L,NKtJE) ) 73.0 
NK=NK+#1 

NK=0 

CONTINUE 

NA=0 

NK=0 

DOAS522. 1=34 2 

DO 523 J=1,KKK 

NA=NA+#1 

BB (NA) =(YYY (I1,NK#1)+VYY(T,NK4+JB) +YYY (1,NK+JE) ) /3.0 
NK=NK+1 

NK=0 

CONTINUE 

NA=0 

NK=0 

bea r9°°C235, 235, 236723 7) p Ree 

JI=8 

JB=9 

JE=17 

GOr 70 534 

JT=6 

JB=7 

JF=13 

GO TO 534 

JI=4 

JR=5 

JE=9 

DOn S322 f=14 0] 

NA=NA+1 

AM (NA) =(M(NK+1) #M(NK4#JB) +M (NK+JE))/3.0 
AW (NA) = (W(NK41) #W (NK4JB) + (NK#J8) ) /3.0 
NK=NK+1 

NB=9 

GOato (235,239, 240,241), KKK 
DO 962 NA=1,3 

WWOCT (NA) =W (NB+2) -W(NB+1) 
WWNON (NA) =W (NB+4) -W (NB43) 
WWDEC (NA) =W (NB+6) —W(NB+5) 
WWHEN (NA) =W (NB+8) -W (NB+7) 
NB=NB+8 

GOaTO 242 

D6"243. NA=1,3 

WWOCT (NA) =W (NB+2) -W (NB+1) 
WWNON (NA) =W (NB+4) -W (NB+3) 
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WWDEC (NA) =W (NB+6) —W (NB+5) 

NB=NB+6 

GOP TO" 28 2 

DO 246 Na=1,3 

WWOCT (NA) =W (NB+2) —-W (NB+1) 

WWNON (NA) =8 (NB+4) -W (NB+3) 

NB=NB+4 

GO TO (248, 248,249,250) , KKK 

WONDH (4) = (WWHEN(1) +WWHEN (2) *WWHEN (3)) 73.0 

WONDH (3) = (WWDEC (1) +WWDEC (2) +WWDEC (3) ) 73.0 

WONDH (2) = (WWNGN (1) *WWNON (2) *WHYNON (3) ) 73.0 

WONDH (1) =(WWOCT (1) +#WOCT (2) + WHOCT (3) ) /3.0 

NB=0 

NK=0 

BORO 254 275% $292) 253) OK KK 

DO 969 NA=1,3 

RRON (NA) =2.0* ( (RET (NK4#2) -RET (NK#1)) / (WWOCT (NB+1) *WWNON (NB#1)) ) 
RPND (NA) =2. 0* ( (RET (NK+3) -RET (NK+2)) / (WWNON (NE+1) +WWDEC (NB+1)) ) 
ERDH (NA) =2.0* ( (RET (NK4+4) -RET {NK+3))/ (WWDEC (NB+1) #*WWHEN (NB+ 1) ) ) 
NK=NK#+4 

NP=NB+1 

GOVTD +254 

pO 255 NA=1,3 

RRCN {NA)=2.0* ( (RET (NK+2) ~RET (NK+1)) / (WWOCT (NB+1) +WWNON (NB+1)) ) 
ERND (NA) =2.0* ( (RET (NK+3) -PET(NK+2)) / (WWNCN (NB#1) +WWDEC (NB+1) ) ) 
NK=NK+3 

NB=NB+1 

GO "TO #254 

DC’ 257 *Na= 1, 3 

RRON (NA) =2.0* ( (RET (NK#2) -RET (NK+1)) / (WWOCT (NB+1) +WWNON (NB+1) ) ) 
NK=NK+2 

NR=NB+1 

GOPTO." (258 7258 259,260)", BKK 

RDH=(RRDH(1) #RRDH (2) +RRDH (3) ) 73.0 
RND=(RRND(1) ¢REND(2) +RRND(3)) 73.0 

FON= (RRON (1) +RRON (2) +RRON (3) ) 73.0 

NB=0 

NK=0 

GO "TO ™ (26°77 26,6202 £203)", Ks 

DO 975 NA=1,3 

THPOCT (NA) =16.0* {RET (NB+1) /WWOCT (NK#1) ) **2 

THPNON (NA) =16.0* (RET(NB+2) /WHNON (NK+ 1) ) **2 

THEDEC (NA) =16.0* (RET (NB+3) /WWDEC (NK4#1) ) **2 

THPHEN (NA) =16.0* (RET (NB+4) /WWHEN (NK+1) ) *¥*2 

NB=NB+4 

NK=NK+1 

GO TO 264 

DO 265 Na=1,3 

THPOCT (NA) =16.0* (RET (NB+1) /WWOCT (NK41) ) **2 

THENON (NA) =16. 0* (RET (NB+2) /WWNON (NK#1) ) **2 

THPDEC (NA) =16. 0* (RET (NB+3) /WWDEC (NK#1) ) **2 

NB=NB+ 3 

NK=NK+1 

GO TO 264 

DO 268 NA=1,3 
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235 


THPOCT (NA) =16.0* (RET (NB+1) /WWOCT (NK+1) ) **2 
THENON (NA) = 16. O* (RET (NB+2) /WWNON (NK+1) ) **2 
NB=NB+2 

NK=NK+1 

GO Woe Wee, 20707271, 272) KKK 

TPCNDH (4) = (THPHEN (1) #THPHEN (2) #THPHEN (3 
TPONDH (3) = (THPDEC (1) + THPDEC (2) #THPDEC (3 
TPONDH (2) = (THPNON (1) +THPNON (2) +THPNON (3 
TPONDH (1) = (THPOCT (1) #THPOCT (2) +THPOCT (3 
GORE RED 7279) 0 702 275 )e RE 

ARONDH (4)=TRE(4) -A 

ARONDH (3) =TRE(3)-A 

ARONDH (2) =TRE(2) -A 

ARONDH (1) =TRE(1)-A 

GOero 4276, 276,277, 278)., KKK 

CPONDIH (4) = (IRE (4) -A) JA 

CPONDH (3)= (TRE(3)—A) /A 
CPONDH (2) = (TRE(2)-A) /A 

CPONDH (1)=(TRE(1)-A) ZA 

NR=0 

NK=0 

COWIOM (281,281, 2625283) KKK 

DO 279 NA=1,3 

SKCCT (NA) = (W(NB+2) -RET(NK#1) ) / (RET (NK+1)—W (NB+1) ) 
SK NON (NA) = (W (NB#4) -RET (NK#2) ) / (RET (NK+ 2) -W (NB#3) ) 
SKDEC (NA) = (W (NB+6) -RET (NK#3) ) / (RET (NK+ 3) —W (NB+5) ) 
SKHEN (NA) = (W (NB+8) -RET(NK#44) ) 7 (RET (NK+4) -W (NB+7) ) 
NE=NBt8 

NK=NK+4 

GO TO 286 

DO 499 NA=1,3 

SKCCT (NA) = (W (NB+2) ~RET (NK4#1) ) / (RET(NK# 1) —W (NB#1) ) 
SKNON (NA) = (W (NB+4) -RET (NK+2) ) / (RET (NK+ 2) -W (NB#3) ) 
SKDEC (NA) = (W (NB+6) ~RET (NK43) ) / (RET (NK+3) -W (NB#5) ) 
NB=NB+6 

NK=NK+3 

GO TO 286 

DO 289 NA=1,3 

SKCCT (NA) = (W(NB+2) ~RET (NK#1) ) / (RET (NK#+ 1) -W (NB+1) ) 
SKNON (NA) = (W(NB#4) -RET (NK4#2) ) 7 (RET (NK+ 2) —W (NB+3) ) 
NE=NB+4 

NK=NK+2 

GO .TO (291,291, 2927293) 7 KKK 

SKONDH (4) = (SKHEN(1) +SKHEN (2) +SKHEN(3)) /3.0 

SKONDII (3) = (SKDEC (1) #SKDEC (2) +SKDEC (3) ) /3-0 

SKONDH (2) = (SKNON(1) #SKNON (2) #SKNON(3)) /3.0 

SKONDH (1) = (SKOCT (1) +SKOCT (2) +SKOCT (3)) /3.0 

GO P10) (297 7.257:7298,295) KKK 

EPONDH (4) =16. 0* ( (TRE (4) -A) /WONDH (4) ) **2 

EPONDH (3)=16.0* ( (TRE (3) -A) /WONDH (3) ) **2 

EPONDH (2) =16. 0* ( (TRE (2) -A) /WONDH (2)) **2 

EPONDH (1) =1€.0* ( (TRE (1) -A) /WONDH (1) ) **2 

GC TO (294,294,295,296) , KKK 

H EONDH (4) = (CLVEPONDH (4) ) 1000.0 

HEONDH (3) = (CL/EPONDH (3) ) *1000.0 
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294 HEONDH (2) =(CL/YEPONDH (2) ) *1000.0 
HEONDH (1) = (CL/EPGONDH (1) ) *1000.0 
NB=0 
GOT LOa(5 05%, 203,305,305) , ARS 
305 DO 306 NA=1,4 
AVGE (NA) = (AREA (NB4+1) +AREA (NB+5) +AREA (NBt9) ) /3.0 
306 NB=NB+1 
GOr TO 307 
304 DO 308 NA=1,3 
AVGE (NA) = (AREA (NB+1) +AREA (NB+4) +AREA (NB+7) ) /3.0 
308 NB=NB+1 
GO TQ 307 
303 POP 309: NAH1,2 
AVGE (NA) = (AREA(NB#1) +AREA (NB+3) +AREA (NB+5) ) /3.0 
309 NB=NB+1 
SOT WET PEG S903) (CARDYT) ,1= 153) 5 (CARD (1) £1=4,6), (CARD (1) ,I1=7,9) 5 
CHOAP DN 146 127042) 
903 FORMAT ('-',133,3A4,T48,3A4,763,3A4,T80, 304) 
NP=0 
DO 904 NA=1,3 
WRKITE(6,905) NA, (AREA (NB4*1) , IT=1,KKK) 
904 NB=NB+KKK 
WRITE (6,906) (AVGE(I) , I1=1,KKK) 
SOS FPORMAT(T l2, SAR BA CCTOTAL, COUNTS) = * 7 
CTVD ,PUNI ECT IGN oy POS, bie L4G yl 1060, LOS, 10.0, T80,F 1050) 
I0G™ POSEAT (T19 ,* AVERAGE AREAS", 
CiPs676r10.0,746, FG. 0, 763781000, 13.0 ,F 10.0) 
NB=0 
DO 920 NA=1,3 
WRITE (6,318) NA, (RET (NB+tI) , 1=1, KKK) 
920 NB=NB+KKK 
WRITE (6,922) (TRE(I) ,1=1,KKK) 
3718 FORMATITI,* RETENTION’ TIME YSECONDS)=*, 
ETSOPSINIECTION PIBVT SES ERIE Ta SPE 7S 24 T64,F7.2,7T81,F 7-2) 
922 FORRAT(T9S; *AVERAGE RETENTION © TINE="; 
CUM GSE] 2 TUS Bley Cuptwas pLOl estas Zc) 
NE=0 
DO 923 NA=1,3 
WRITE (6,924) NA, (PEAK (NB4+1) , 1=1,KKK) 
923 NB=NB¢tKKK 
WRITE (6,925) (KPEA(I) ,I1=1,KKK) 
S2U FORMAT (T11, . PEAK HEIGHT COUNTS) =*, 
CE59 SAN IECTION!;, 0370 0808 ne0 yt ey 0, Te4, £720,161, Fi .0) 
O25) FGRNAT(T 12, tAVERAGE PEAK HEITGHT=", 
CIBSEFE).O TUG Fi. O, TORE OT oie ry. G ) 
NB=0 
DO 926 NA=1,3 
WRITE (6,927) NA, (XXX (1,NBtL) ,I=1, KKK) 
926 NB=NB+KKK 
WRITE{6,928) (AA(I),1=1,KKk) 
927 FORMAT (TIE, XT "COORDINATE =*7 
CTOTPUTINIECTION PL sel 14, Poe eo te 25T64, FI. 2,78 1,6 1< 2) 
928 FORMAT (T10,* AVERAGE X7*°CCORDINATE='; 
CTBST FI «2, TESTE la 2g LON pt eg Ol ghilee) 
NB=0 
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GC TO{34C0,340, 341, 342) , KKK 
JF=8 
GOnTOnaleAv 
341 JB=4 
JE=6 
GO TO 1342 
340 JR=3 
JF=4 
1IZS27 E929 NA=1,3 
WRITE (6,930) NA, (XXX (2,NBt1) ,1=1, KKK) 
929 MB=NB+KKK 
WKITE (6,931) (AA (I) ,1=JB,JE) 
930 FORMAT (T18,*X2 COORDINATE=', 
Cl eh COA ONY LO TIN ele epi dey lee yO, Fis eet ly fied) 
Dap GEMAT ({T10,*' AVERAGE X2 COORDINATE=", 
Cl ae oe LODE le 25 Ob ti ae yo ler lsc) 
NR=0 
GO TO (348,348, 349, 350) , KKK 
350 JB=9 
JF=12 
GOTO: 113/50 
349 JR=7 
JE=9 
GO TO 1350 
348 JB=5 
JR=6 
T1350) 2007 932 NA=1,3 
WPITE (6,933) NA, (XXX (3,NB4I) ,I=1, KKK) 
932 NBE=NB+t+KKK 
WRITE (6,934) (AA (I) , 1=JB,JF) 
OSseFORGAT (T138,"°X3 COGRDINATE=", 
CPOs, “NI EGON! (23 ,tSe Tie, LOO, Fie 2, lt OU, Fl. 2,18 1,F 7.2) 
934 FORMAT (T10,"AVERAGE X3 COCRDINATE=', 
OT 34 Bieciguens ph le Ze 104, F Pe2, TS 1 gt lee) 
NR=0 
GOTO: , (3/5673 36935 17,7358 eK KR 
358 JB=13 
JE=16 
GO M0 w1p.58 
357 JB=10 
JE=12 
GOO; 1358 
356 JB=/7 
- JF=8 
1358 DO $35 NA=1,3 
WRITF (6, 93€) NA, (XXX (4U,NBtI) ,I=1, KKK) 
935 NE=NB+KKK 
WRITE(6,937) (AA (I) , I=JB,JE) 
936 FORMAT. (T18,'X4 COORDINATE=', 
CT 9S. eiN ECT iG Nass T See fee GF 7.2 TOU, Fs Jel 81, F le 2) 
937 FORMAT(T10,"AVERAGE X4 COORDINATE=!, 
CRI Tief dig EUG Eilat flO Wi) 22), TS Waris 62) 
NB=0 
DO 938 NA=1,3 
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WRITE (6,939) NA, (YYY(1,NB+I) ,1=1,KKK) 
938 NB=NRtKKK 
WRITE (6,940) (BB(I) ,I=1, KKK) 
939 abOBMAT(T 11, °Y 7 AND vu COORDINATE=', 
CT99,"INJECTION',13,133,F8.1,T48,F8.1,T763,F8. 1,780, F8. 1) 
940 FORMAT (T3,'"AVERAGE Y1 AND Y4 COORDINATE=!", 
CT33,F8.1,T48,F8.1,7T63,F 8. 1,780,F 8. 1) 
NRE=0 
BOT Oa Iine 3h2 yd uit TH)» KEE 
374 Jb=5 
JE=8 
GO. TO, 1374 
373 JB=4 
JE=6 
GO TO 1374 
3a 3 B=3 
JE=4 
L3vH DO 941 NA=1", 3 
WRITE (6,942) NA, (YYY (2,NB+I) ,I=1,KKK) 
941 NB=NBt+KKK 
WRITE (6,943) (BB (I) , I=JB,JE) 
SUZ. HORMATATI1,°Y2 AND Y3 COORDINATE=', 
CTP sp LNJECTION?, 13,%33,F8.1, 748, F8.41,163,F8.1,7T80,F8.1) 
943 FORMAT (T3,'AVERAGE Y2 AND Y3 COORDINATE=", 
sg bs Vp tte, Ue dg LOD ,F Os 1p TOO ye oe 1) 
NP=0 
GO) TR. (3.559.396 Say BOB) ok KK 
398 JB=7 
GO. TO. 13,98 
3974 SBEH5S 
GOTO 31.398 
396 JP=3 
1398 DO 950. NA=1,3 
WRITE (6,951) NA, (M(NB+I) ,I=1,JB,2) 
950 NE=NBtJB+1 
WRITE (6,952) (AM(I) ,T=1, JB, 2) 
IOAN FORMAT (Til? POSITIVE, SLOPE=3 , 
CT99,*INJFECTION', 13,734, F7.0,T49, F7.0,T64,F7.0,T81,F7.0) 
952 FORMAT {(T9,*AVERAGE POSITIVE ILOPE=*, 
CT34,F7.0,T49,F7.0,T64,F7.0,781,F7.0) 
NB=0 
GO TO (404, 404,4C5,406) ,KKK 
406 JE=8 
GO TO 1406 
405 JB=6 
GC TO 1406 
4O4W JB=4 
1406 DO 953 NA=1,3 
WRITE (6,954) NA, (M(NB+I) ,I=2,JB,2) 
953 NB=NB+JB 
WRITE (6,955) (AM (I) ,1=2,JB, 2) 
954 FORMAT (T17,'NEGATIVE SLCPE=/., 
CTO of INJECT LONI dy te hy Fs oO, THD, Fe) Dy TON, Fi) 00% TS lg Pal 0) 
955 FORMAT (T9,*AVERAGE NEGATIVE SLCPE=', 
CT34,F7.0,T49,F7.0,T64,F7.0,781,F720) 
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NB=0 
GO TO (412,412,413,414) , KKK 
414 JB=7 
GO TO 1414 
413 JB=5 
cor To 1474 
412 JB=3 
1414 DO 956 NA=1,3 
WRITE (6,957) NA, (W(NB+I) ,T=1,JB,2) 
956 NR=NB+JB+1 
WRITE (6, 958) (AW( 1) 121708, 2) 
957 FORMAT (T4,"POSITIVE BASELINE INTERCEPT=', 
CTEOG Ss ENISOCINON® ,13,738, 77.2, 795, F7,2,T64, 702,181, 87.2) 
958 FORMAT (T5,'AVERAGE tVE BASE INTERCEPT=', 
CMa OF FTS 240495 £9, 22 6ue PP .2,781,F 7.2) 
N B=0 
GO! TO? (4:20, 4205421, 422)), KKK 
422° OB=8 
Ge 10'1422 
421 JB=6 
GOMmesaho? 
420 JB=4 
1422 DO 959 NA=1,3 
WRITE (6,960) NA, (W(NB#T) , T=2, JB, 2) 
959 NB=NB+JB 
WRITE (6,961) (AW(I) , 1=2,JB, 2) 
960 FORMAT (T4,*NEGATIVE BASFLINE INTERCEPT=", 
CTO9, INJECTION! 13,734, 97.2,1T49, 67. 2,164, F7.2;781,F7<2) 
961 FORMAT(T5,'AVERAGE -VE BASE INTERCEPT=', 
CTRL, FT. 2), TURF .2LTSRe FC ol, TENA FASO) 
NB=0 
GO TO (428,428,429,430) ,KKK 
430 DO 966 NA=1,3 
WRITE (6,967) NA,WWOCT(NB+1) ,WWNON (NB+1) ,WWDEC (NB+1) ,WWHEN (NB#1) 
966 NB=NB+1 
GO TO 431 
967 FORMAT (T12,*BASELINE PEAK WIDTH=", 
CT99,? INUVECTHON?), 13% TRAG Pe 24k Foo, 764, F7. 2,781, F 7s 2) 
429 DO 432 NA=1,3 
WRITE (6,967) NA,WWOCT(NB+1) , WWNCN (NB+1) , WHDEC (NB¥1) 
432 NB=NB+1 
GOA TOR BOH 
428 DO 434 NA=1,3 
WRITE (6,967) NA, WWOCT(NB+1) , WWNON (NB+1) 
434 NR=NB+1 
431 WRITE(6,968) {WONDH (I) ,I=1,KKK) 
968 FOPMAT(T4,*AVERAGE BASFLINE FEAK WIDTH=", 
Cnet U2, 149,577. 2,164 rie, Th 1 FT. 2) 
NB=0 
GC TO (436,436,437,438) , KKK 
438 PO 972 NA=1,23 
WRITE (6,973) KRON(NE*1) , RRND(NB#1) ,RRDH (NB+1) ,NA 
973) FORMAT (T21, *RESOLUTION=! ,TH2,F7.2,,1T57,F7c2,T72yF7.2,7T99, 
C'INJECTION',13) 
972 NB=NB+1 
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WRITE (6,974) RON, RND, RDH 
774 FORMAT (T13,"AVERAGE RESOLUTION=! ,T42,F7.2,T57,F7.2,T12¢F7«2) 
GO TO 439 
“37 DO “io NAST 3 
WRITE (6,441) RRON(NE+1) , RRND(NB+1) ,NA 
441 FORMAT (T21,*RFSOLUTION=",T42,F7.2,T57,F7.2,T799," INJECTION',13) 
4W4UC NB=NB+1 
WRITE (6,974) RON, RND 
GO TO 439 
436 DO 442 NA=1,3 
WRITE (6,443) RRON(NE+1),NA 
443 FORMAT (T21,"RESOLUTION=',142,F7.2,T99,*INJECTION',I3) 
442 NB=NB+1 
WRITE(6,974) RON 
439 NB=0 
GO TO (444, Guu, 445,446) , KKK 
446 DO 979 NA=1,3 
WRITE(6, 980) NA, THPOCT (NB+1) , THENON(NB+1) , THPDEC (NB+1) , THPHEN (NB+1) 
979 NR=NB+1 
GO TO 407 
980 FORMAT (T13,"THEORETICAL PLATES=!, 
C1927, LNJECTION* 13,034, F7.0,149,F7.0,164,F7.0,T81, FT. 0) 
445 DO 448 NA=1,3 
WRITE (6,980) NA, THPOCT (NB+1) , THPNON (NB+1) , THP DEC (NB#1) 
448 NB=NB+1 
GO TO 447 
444 DO 450 NA=1,3 
WRITE (6,980) NA, THPOCT (NB+1) , THPNON (NB+1) 
450 NB=NB+1 
447 WRITE (6,981) (TPONDH(I) ,1=1,KKK) 
981 FORMAT(T5,"AVERAGE THEORETICAL PLATES=', 
CT34,F7.0,T49,F7.0,1T64,F7.0,781,F7~0) 
WRITE(6,982)A 
982 FORMAT(T21,'GAS HOLDUP=',F7.2,2X,"SECONDS') 
WRITE (6,983) {ARONDH(I) , I=1,KKK) 
983 FORMAT(T8,'ADJUSTED RETENTION TIME=", 
CEI AVERAGES 230, 0962) [toga fou, rye 2, Tol, E72) 
WRITE (6,984) (CPONDH (I) ,1=1,KKkK) 
984 FORMAT (T4,*CAPACITY OR PARTITION RATIO=!, 
CIO? 7 AVERACE, T34 Fi. 2, 0h Vesa oe tod ers. 2,18 1,87 .2) 
WRITE (6,467) (HECNDH (I) ,1=1,KKK) 
467 FORMAT(T21,*HETP (MM) =!?, 
C199, "AVERAGE, T34. ,F7.2,149,F7.2,764,,F 7. 2,78 1, F 7-2) 
WRITE (6,474) (EPONDH (I) , 1=1,KKK) 
474 FORMAT (T9,*EFFECTIVE PLATE NUMBER=', 
C99, AVERAGE! ,T397 Fo 507499 092.0 164, £7 60,1815 £720) 
NP=0 
GO TO (478&,478,479,480) ,KKK 
480 DO 481 NA=1,3 
WRITE (6,482) NA,SKOCT (NB#1) ,SKNON (NB+1) , SKDEC (NBt 1) ,SKHEN (NB#1) 
481 NB=NB+1 
GC TO 1480 
482 FORMAT (T27,"SKEW=", 
CT97, NIECTION* 13,2357 50.4 yo, F7.4 | T64, PI. 4,781, FT) 
479 DO 485 NA=1,3 
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WRITE (6,482) NA, SKOCT(NB+1) ,SKNCN (NB+1) , SKDEC (NB#1) 
485 NB=NBe1 

GO TO 1480 
478 DO 487 NA=1,3 

WRITE (6,482) NA,SKOCT (NB+1) , SKNON (NB+ 1) 
487 NB=NB+1 
1480 WRITE (6,483) (SKCNDH(I) ,1=1,KKK) 
483 FORMAT(T19,"AVERAGE SKEW='", 

CTP348,F7.4,T49,F7.4,764,F7.4,1781,F7.4) 

900 CONTINUE 

STOP 

END 
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REAL T,M1,"2,TC1,TC2,VC1,VC2,DGO0,DL,L,K,ALPH,J,TK1,TR2,EXPR1,H, 
CPERC,S,ROHS, PL,PO,FR, PR,X1,X2,X3,X4,Y¥1,Y2,Y¥3,Y4, EXPR2, THEORR 
Coy 4 

INTEGER JI 

READ (5,700) S,ROHS,M1,M2,TC1,TC2, VC1,VC2 

70Q FORMAT (8F10, 3) 

DOP GOO. T=4 100, 4 

READ (5,910) TR1,TR2,FR,PI,PO,T, PR 

GAC SEORMAT (7 F70.3) 

WRITE(6,911) 

911 FORMAT (1 0", 29,°5", 117 ,.) ROHS" 128, 'TR1*,7T37,1TR2* P47, ORR? 
CT98,*PI*,T68,'PO', 78,9 PR?) 

WRITE (6,912)S,ROHS,TR1,TR2,FR, PI, PO, PE 

pe eT OEMAT (4708 8F 10.3) 

WRITE(6,921) 

924 POEBAT(V OM, TL, * MA pad eM 2 9268 TC1 tyr I6 MT CAA TaG, ty CTs, 
Cro, Vices) 

WRITE (6, 922) 01582 TC 1 Teo mvc Ay vc? 

922 FORMAT (6F10.3) 

ALPH=TR2/TR1 

J=1.5* ((PI/PO) **2-1.0) 7 ( (PI/PO) **3-1.0) 

R1=0.. 434 (T4273. 0) 7100.0) **4 28 1 

A2=((1.0/781)4(1.0/M2) ) **0.5 

X3=0.92* ( (TC 1*TC2) 710000. 0) **0.1405 

X4=(((VC1/100. 0) ** 0.4) + ((VC27/100.0) #0. 4) ) #*2 

DGO= (X1*X2) / (X3*X4) 

DL=0.00001*DGO 

READ (5,81C) L,PERC,K,EXPR1,H 

SAO “FORMAT (5510.3) 

Y1=0.25*K* ( (ALPH-1.0) /ALPH) *L**0.5 

Y¥2=2*J *0.6* TR2*DGO* (1.0+#K) /L 

Y¥3=0. 667¥*L¥* (1.04+K) **0.5*PERCKK2 

Y4=DL*TR2*K* (S*¥ROHS) **2 

Z1=(TR2+TR2**0. 9) /TR2 

YS=ALOG10 (1. 0+#L/15) * (2®K**Z1) /K¥*0. 7 

THEORR=Y 1/ (Y24+Y3/Y44+Y5) **0.5 

EXPR2=0.25* ( (ALPH-1.0) /ALPH) * (K/ (K#1.0) ) * (L/H) **0.5 

WEITT:(6, 93.0) 

930 FORMAT ('0",T6,"ITEMP',T18,'DGO',T28,'DL* ,T38,'J*,T47,*ALPH', 
CT58,'K"',T64,"PER CENT',T76,"COL L*,T85,!THEOR R',1T97,*EXP R1*, 
ClIUG, Bsr R274 

WPITE (6,811) T,DGO,DIL,J,ALPH,K, PERC,L,THEORR, EXPR1,EXPR2 

CUS FORMAT (* (*, F9.37P905,2%,F 10.9, 6F 10.3} 

900 CONTINUE 

STOP 

END 


phic: © HAA ai i 
¥,S9 PMs ee: 


, ~ : i | 
SOW POY SOT, FIT, SH PR 
i 7 » 


‘Dees 


ATT OT.TI 


| ss" 
rat’. ost, '2non* TPP AP ae 
cay fer tamed 


4 lia wae ; 
3 ee i: {' 
ae - he 
; ~~ The 
= 
; 


i 


1O¥) 
¢ 


? 


TPT RMTs reT, e8O 5.8 4 re | y) ™ 


, Si 4% re *0. 7 


€ 01% 


_ 


“a 
5 0 a +h 


oa | 4 ar 


Vf Vv, SOT, FOR. nin esse ayaTIT 


re oe 
rar\ 


My\{9 t-sestonta) seats 


. 

f re (0, | 
# 

PEs iw 


0. 1) aGONNSAT HA. Or Les 


Ory (0 ETS +R) * 

eLESUMN OLE) 4 + ePeyO.e )* 

OGOTS OTM SEY) @ .92€ 
‘(7.998 (06 OOP FOV spa 

aXxVE XS ESKFTA 000 

oager9900,0 

TY 3. Hs i%, 3 (OTs, ae 

fonda) THR 

, ti nias} ound 0 


bo0s 


wt 


9 OFS INEM s ry eevee 088 y 


c+#* (2HO8* ay ¥Heg aT? 
AT (2 .O** SATS AT yar 
é 1 * (Bt S+0 of) OF Bom 
*# (PY oR ¥XEY SS ¥p\F t=2 - 
nt BI EADE SES a GeSs 
(O€2 a) OTF 
rT, "Onatr? .et 03) DANA 
tra tat  baRs 40 e8 dy 
(‘<a 43a" Riite 
L, 10, 000. TLPRR, 4) YTT 9 
By Pe PG yes eq, ‘yA iOF 
= ire bb 
To} "¢ 


aif 


— a 


Bl 
* 
: \ ) 
£ } 
P \ 

i olan ay } 

Ue lah 

ee 
-) - y 

‘} 

t. \ 


; ere 
72) 


i, >: if 
ie. Ly lg 
a i 
May aan 


Qs es eh onmeents @: 
Sin iee gate een Ste oN 
peered 


ap nee 
sevan? ioe jase wae 
¢ sp Aiold 


Lea a 


ys 


iat tel 


ees 


Reis ihea aly 


Les eoe. 
oh errno 

mb edhe sasetinn 
vad gascaapate 


pep oea re) 


rar Pere ies eee ene STs 


Merete re tere tear acer a 


Lotiows Lat mn amelie 
apeimye ec sh) em 


